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Abstract 

A method of optimisation of new extractants structure using the desirable function has been developed. Earlier the 
desirable function has been proposed by Harrington (Ind Qual Control 21: 494-498, 1965) for the optimisation of 
processes with several response functions. The developed method of optimisation of new extractants structure has 
been used for construction of phenolic type extractants (PTE) (a class of N-(2-hydroxy-5-nonylbenzil)-dialkylamines). 
It has been offered to use the charge on the nitrogen atom, the heat of dissociation of phenolic group, the 
logarithm of distribution factor of extractant between water and octanol (computed data) and maximum 
permissible concentration of extractants in aqueous phase (MPC) of the o-replaced phenols (the literary data) as the 
controllable parameters, defining efficiency of extractants for molecular design of PTE. During optimisation of 
extractants structure the quantity of alkyl substitutes at nitrogen atom, the carbon atoms number in these 
substitutes and the electronegative substitutes in o-position to phenolic group have been varied. As the result of 
the molecular design, the optimal structure of PTE found is N-(2,3-dihydroxy-5-nonylbenzil)-didecylamine, which 
perfectly meets the requirements to industrial extractants. 
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Introduction 

The search and development of new effective extractants 
are an important scientific and practical issue (Zolotov 
1982). There are some requirements to the extractant 
for metal extraction when using the solvent extraction 
method. However in practice, it is difficult to find an 
extractant responsible for all requirements, so one usu- 
ally compromises and chooses a solution (Ritcey & 
Ashbrook 1979). The testing of substances, used in other 
industries as extractants, is of big importance in the ex- 
pansion of the range of extractants for the extraction of 
metals (Zolotov 1982). We suggest the use alkyl deriva- 
tives of phenols, applied in synthetic rubber industry, 
plastic industry and elastomers industry, as extractants 
(Bukin et al. 1999; Bychenkov et al. 2008; Gladikova 
et al. 2002). However, despite the promise of this idea, it 
is necessary to consider the fact, that reagents developed 
for other purposes may not meet to all requirements for 
extractants. That is why the problem of extractant mol- 
ecule design, that meets the certain requirements, arises. 
The computer technology development in recent times al- 
lows solving this problem by using calculation methods, 
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like the methods of quantum chemistry or the method of 
group contributions, etc. The development of a new, more 
effective extractant allows to reduce expenditures for the 
rare metals extraction and opens new perspectives in the 
certain metals usage, that are in limited use because of 
their high price (for example, scandium) (Korshunov et al. 
1987; Komissarova 2006). 

In the works (Hay 2008; Varnek 2008) in the course of 
the molecular design of extractants thermodynamic cri- 
teria are being considered: the energy of interaction be- 
tween the receptor-ion (Hay 2008) and the constants of 
solvent extraction, distribution coefficients and separ- 
ation coefficients of extractable metals (Varnek 2008), 
with the non thermodynamic factors such as MPC 
(maximum permissible concentration) are not taken into 
account. 

The purpose of the current investigation is the develop- 
ment of a method for optimisation the structure of new 
extractants using the desirability function, proposed earlier 
by Harrington (1965) for the optimisation of processes, 
characterised by a few response functions (Ahnazarova & 
Kafarov 1978). 

According to Harrington, the desirability function is a 
dimensionless scale, which allows to convert any 
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Figure 1 N-(2-hydroxy-5-nonylbenzil)-dialkylamine (NBAA); R : 

Qi^im-I- 



response so that it is interpreted in the terms of useful- 
ness or desirability for any specific application. 

For the unilateral restrictions like y < y max or y > y min 
(y-response function) a suitable form of the transform- 
ation y in d (particular desirability function) is the expo- 
nential function: 



d = exp[—exp(—y')} 

where: 

y = b 0 + hy 



(1) 



(2) 



The coefficients b 0 and bj will be determined, if for 
the two values of given property y one sets the corre- 
sponding values of the desirability d, preferably in the 
range 0.2 < d < 0.8. 

Having a few responses, converted into d scale, 
it is possible to combine the generalised desirability 



Table 1 The values of q and logP for NBAA extractants at 
varying lengths of two hydrocarbon radicals 

Extractant 

NBAA-00, n = 0 
NBAA-11, n = 1 
NBAA-22, n = 2 
NBAA-33, n = 3 
NBAA-44, n = 4 
NBAA-55, n = 5 
NBAA-77, n = 7 
NBAA-1010, n = 10 



OH 




N 



' C 10 H 21 
^ C 10 H 21 



C 9 H 19 

Figure 2 The substituted in the o-position N-(2-hydroxy-5- 
nonylbenzil)-didecylamine. 



index D from the different d, using the following 
expression: 



D = \J d\di ■ ■ ■ d~k 



(3) 



Methodology 

The proposed method of the optimisation the structure 
of new extractants has been used for the designing 
phenolic type extractants (PTE) (class N-(2-hydroxy-5- 
nonylbenzil)-dialkylamines). Moreover the following 
controllable parameters were selected: 

1. The formal charge on the nitrogen atom (q), as 
some metals such as scandium (Bychenkov et al. 
2008; Gladikova et al. 2002) are extracted by PTE 
with the formation of chelates and the nitrogen 
atom inclusion to the chelate cycle. The value of 

q was calculated by the Milliken approach with the 
DFT (B3LYP) method, 6-21G basis set using the 
GAMESS-US (version September 7, 2006) program 
package. 

2. The logarithm of the partition coefficient for 
n-octanol/water (logP), calculated by means of the 
group contribution method using the program 
Chem3D Ultra version 7.0.0, included in the package 
of applied programs (PAP) ChemOffice Ultra 7.0.1. 
The value logP allows estimating the solubility of the 
extractant into the aqua phase and, therefore, the 



0.093 


logP 

5.11 


Table 2 The dissociation energy AE D , 
nitrogen atom and logP of NBDA- R' 


the charge on the 
extractants 


0403 


5.99 


Extractant 


A£, a.u. 


<7 


logP 


0449 


7.05 


NBDA - OH 


0.598 


-0.285 


14.9 


0443 


8.11 


NBDA - N0 2 


0.551 


-0.279 


13.8 


0440 


9.16 


NBDA - F 


0.582 


-0.278 


15.8 


0436 


10.22 


NBDA - CI 


0.587 


-0.294 


16.2 


0428 


12.34 


NBDA - Br 


0.583 


-0.294 


16.4 


0.456 


15.51 


NBDA - 


0.576 


-0.279 


16.6 
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Table 3 The values of o-substituted phenols MPC 


Substitutes in ortho-position to 


MPC W (in aqua phase). 


phenolic group 


mg/m 3 


without R' (Chernyshev et al. 2005) 


0.001 


-CI (Chernyshev et al. 1 999) 


0.0001 


-N0 2 (Chernyshev et al. 2004) 


0.06 


-OH (Chernyshev et al. 2004) 


0.1 



possible losses of the extractant during the solvent 
extraction process. 

3. The energy of dissociation of the phenolic group 
{AE d ), calculated by the DFT (B3LYP) method, 6- 
21G basis set using the GAMESS-US (version 
September 7, 2006) program package. According to 
(Gladikova et al. 2002), in the course of scandium 
solvent extraction by the N-(2-hydroxy-5- 
nonylbenzil)-(3,p-dihydroxyethylamine (NBEA) 
extractant the optimum pH is 4.5. However, in 
industrial products containing scandium there is a 
large amount of iron (III) as a rule (Korshunov et al. 
1987), and during the precipitation of iron 
hydroxide large losses of scandium take place due to 
co precipitation. So, the problem of PTE designing, 
which extracts scandium at lower pH, arises. 
Therefore, the choice of given option is specified by 
the necessity of increasing the acidity of the 
phenolic group by introducing electronegative 
substituents in the ortho-position to it (Nesmeyanov 
& Nesmeyanov 1970). 

4. MPC of o-substituted phenols (Chernyshev et al. 
1999; Chernyshev et al. 2004; Chernyshev et al. 
2005). 

As an initial structure for further optimisation N-(2- 
hydroxy-5-nonylbenzil)-dialkylamine (NBAA) was selected 
(Figure 1). 

Earlier, the representatives of this class of 
extractants - NBEA and N-(2-hydroxy-5-nonylbenzil)- 
|3-hydroxyethylamine (NBEA-1), which are used in the 



synthetic rubber industry, were tested for scandium 
solvent extraction. 

Results and discussion 

At the first stage of investigation with the purpose of 
extractant structure optimisation the length of hydrocar- 
bon substituents R was varied at the nitrogen atom. The 
maximum length of the hydrocarbon radical was limited 
to ten carbon atoms in order to avoid the high viscosity 
of the extractant, the decrease of its capacity of extracted 
component and the increase of the steric hindrances of 
complexation (Rozen & Krupnov 1996). The results are 
shown in Table 1. 

As can be seen in the Table 1, the hydrocarbon radical 
length increasing both the charge of the nitrogen atom 
(in absolute value) and the logarithm of the distribution 
of NBAA between the aqua and octanol rises. The in- 
creasing of the charge at the nitrogen atom leads to rise 
in extraction ability of this extractant, and an increase of 
logP shows a decrease of its solubility in aqueous phase. 
Obtained results are consistent with the data given in 
(Rozen & Krupnov 1996). Consequently, in the first 
stage of optimisation NBAA-1010 was the most effective 
extractant. 

In the second stage of the structure optimisation of 
phenol type extractants by varying of the substituents in 
the ortho-position to the phenolic group o-substituted 
N-(2-hydroxy-5-nonylbenzil)-didecylamine was used as 
the original compound (Figure 2). 

The acid dissociation constant of the phenolic group 
was evaluated by the o-substituted phenols dissociation 
energy, which was calculated by DFT (B3LYP) method, 
6-21G basis set. 

The calculated values of the dissociation energy AE D 
for NBDA-RJ the charge on the nitrogen atom and the 
logarithm of the distribution coefficient between the 
aqua and octanol are shown in Table 2. It is obvious that 
the charge on the nitrogen atom changes to a small ex- 
tent probably due to the remoteness of the electronega- 
tive substituent relatively the nitrogen atom, and the 



Table 4 The calculation of the generalised desirability function 



Substitutions in 
o-position 


Particular desirability functions 






Generalised desirability functions 


d, (q) 


d 2 (logP) 


ds (AF D ) 


d 4 (MPC w ) 


D? 


D*P 


Without the substituent 


0.690 


0.689 


0.400 


0.290 


0.579 


0.476 


-N0 2 


0.430 


0.700 


0.500 


0.500 


0.494 


0.505 


-F 


0.400 


0.726 


0.560 




0.546 




-CI 


0.700 


0.771 


0.530 


0.250 


0.659 


0.495 


-Br 


0.690 


0.792 


0.550 




0.669 




-1 


0.410 


0.810 


0.590 




0.581 




-OH 


0.530 


0.600 


0.460 


0.650 


0.527 


0.541 



* D 3 was calculated by the equation D3 — tfd\ x c/2 x d$. 
** D 4 was calculated by the equation D 4 — -tfd-\ x cfa x fife x cf 4 . 
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logarithm of the distribution index by the transition 
from the nitro-group to the iodide ion increases, which 
is consistent with the decrease of the hydrophilicity of 
these substituents. 

An important indicator, determining the possibility of 
the extractant use in rare metal technology, is the max- 
imum permissible concentration (MPC) of extractant 
in the aqueous phase. Considering the solubility of 
extractant in the aqueous phase and its losses in the for- 
mation of an emulsion, one should always take into ac- 
count MPC of the extractant when selecting. During the 
study of literature only 4 MPC values were found, 
namely: for phenol, o-chlorophenol, o-nitrophenol and 
o-hydroxyphenol (pyrocatechol) (Table 3). 

In solving the PTE optimisation problem as responses 
one used: yi = q - the charge on the nitrogen atom, 
j2 = log/ 5 - logarithm of the PTE distribution index be- 
tween water and octanol; j3 = AE D - the energy of the 
phenol group PTE dissociation, a.u.; = MPC W - the 
maximum permissible concentration in water, mg/m 3 ; 
D- the generalised desirability function. 

For the comparative evaluation of the extractants ef- 
fectiveness with various substituents in o-position to the 
phenolic group the generalised desirability function was 
determined by the formula: 



Da 



\fd\ x ~di x ~d% 



(4) 



where dy d% dz, d4 - particular desirability functions. 

It is necessary to convert obtained PTE characteristics 
into dimensionless uniform scale y, to construct particu- 
lar desirability functions. The developed extractant is to 
meet the specified requirements by the four indicators. 
Based on these requirements the values yj, ya y& y& cor- 
responding to two base benchmarks on the desirability 
scale, were selected. 

The response y into the desirability function d was 
converted by the Equation (1). The results are given in 
Table 4. 

Table 4 shows that the extractant with bromine sub- 
stituent in the o-position to the phenolic group has the 



maximum value of the desirability function by three 
criteria (the charge on the nitrogen atom, logarithm of 
the PTE distribution index between water and octanol; 
the energy of the phenolic group dissociation). N-(2,3- 
dihydroxy-5-nonylbenzil)-didecylamine is the most ef- 
fective extractant by the four criteria (besides the three 
above mentioned, also the MPC value in the aqueous 
phase) (Figure 3). 

Conclusions 

Eventually, as a result of the molecular design of the 
extractant there have been found the optimal structure, 
that is N-(2,3-hydroxy-5-nonylbenzil)-didecylamine, which 
meets to the maximum extent to the considered require- 
ments to the industrial extractants. We plan to introduce 
the proposed extractant into the technology of scandium, 
rhenium, gallium, cobalt and other rare metals extraction 
after the synthesis and experimental verification. 

Nomenclature 

d Particular desirability function 

D generalised desirability index 

logP the logarithm of the partition coefficient for 

n-octanol/water 
q formal charge on the nitrogen atom 

y response function 
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